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Abstract. We investigate possible formation sites of the cannonballs (in the gamma ray bursts context) by calculating their 
physical parameters, such as density, magnetic field and temperature close to the origin. Our results favor scenarios where 
the cannonballs form as instabilities (knots) within magnetized jets from hyperaccreting disks. These instabilities would most 
likely set in beyond the light cylinder where flow velocity with Lorentz factors as high as 2000 can be achieved. Our findings 
challenge the cannonball model of gamma ray bursts if these indeed form inside core-collapse supernovae (SNe) as suggested 
in the literature; unless hyperaccreting disks and the corresponding jets are common occurrences in core-collapse SNe. 
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1. Introduction 

It has been argued in the literature that, as an alterna¬ 
tive to the fireball scenario (e.g. Piran I1999I and references 
therein), the so-called cannonball (CB) model provides a 
good fit to the observed GRB flux and temporal variations 
l |Dar & De Riljula 2004^ . For example, to explain GRBs, CBs 
must be created in supernova explosions and accelerated to 
high Lorentz factors, Fcb ~ 1000. However, the origin of 
these highly relativistic “balls” of matter has not yet been in¬ 
vestigated and is the subject of much debate and controversy. In 
order to shed some light on the still open questions of their for¬ 
mation and early evolution we investigate, in this paper, the CB 
physical conditions at the origin given their features at the dis¬ 
tance when they become transparent to their enclosed radiation 
as required to explain GRBs. Our proposal is that the conditions 
within the CB as we scale the distance down along its path to 
the origin should be an indication of their formation site. This, 
despite the simplicity of our approach, we hope might help elu¬ 
cidate some questions related to the origin/existence of these 
CBs. We start in Sect. 2 by a brief introduction to the CB model 
as described in Dar & De Rujula (2004). In particular CBs con¬ 
ditions at infinity which best fits GRB lightcurves are isolated. 
In Sect. 3 we present the methods we adopted to extrapolate 
back to the CB source. In Sect. 4, given the conditions at the 
origin, we study possible formation sites and explore formation 
mechanisms. Sect. 5 is devoted to the study of mechanisms ca- 
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pable of acceleration CBs to Lorentz factors as high as ~ 1000. 
We summarize our results and conclude in Sect. 6. 

2. The CB model for GRBs 

In the CB model for GRBs, the prompt gamma ray emission 
is assumed to be produced when ambient light from the su¬ 
pernova* is Compton up scattered by the electrons in the CB. 
These CBs move with Lcb ~ 1000 with respect to the su¬ 
pernova remnant and as such the emitted radiation is highly 
beamed in the observer frame (Dar & De Rujula 2004). For 
the first 10'* s in the CB rest frame, the CB in a fast cooling 
phase emits via thermal bremsstrahlung (Dado et al. 2002), 
but eventually, it is argued, its emissivity is dominated by syn¬ 
chrotron emission from ISM electrons that penetrate it. A CB 
will become transparent to the bulk of its enclosed radiation in 
a time of (9(1) s in observer frame after it exits the transparent 
outskirts of the shell of the associated SN. The internal radia¬ 
tion pressure drops abruptly and its transverse expansion rate 
is quenched by collisionless, magnetic-field mediated interac¬ 
tions with the ISM (Dado et al. 2002). 

Typical values for CB parameters as derived by Dado et al. 
(2002) and Dar & De Rujula (2004) are given in Table^where 
we denote the radius of the CB by Rqb, the distance travelled 
by the CB from its origin by D, the time passed in the CB 
rest frame by t. The expansion velocity of the CB is denoted 

’ The wind from the SN progenitor star is ionized and is semi¬ 
transparent to photons in the visible and UV frequencies. 
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Table 1. Cannon ball parameters as given by Dado et al. (2002) 
and Dar & De Rujula (2004) 


parameter 

value 

.^CB,max 

2.2 X lO^"* cm 

.^CB,trans 

10^® cm 

.C^trans 

1.7 X 10^® cm 

Fcb 

1.0 X 10® 


1.0 X 10® 

Z 

1 

/3i 

1/V3 

Ncb 

10®° 

Mcb 

10"® g 


by Pi = Vi!c, the number of baryons and the mass of the CB 
by iVcB and Mcb, respectively. The subscript ’’trans” refers 
to the point where the CB becomes transparent (’’transparency 
radius”). 

By fitting the observed GRB afterglow the CB Lorentz fac¬ 
tor is estimated to vary between Fcb = 250 and Fob = 1600, 
while the number of baryons is of the order of 10®°. With this 
information at hand, our goal is to derive the conditions in the 
CB as we integrate back to a plausible source. As an indication 
of the close proximity to a compact source (e.g. black-holes and 
neutron stars) when applicable we will make use of the notion 
of light cylinder which we take to be about i?L ~ 10^ cm. 

3. Cannonball propagation and evolution 

In this section, we will explore the evolution of the CB by ex¬ 
trapolating backward from the location where the GRB occurs 
to where the CB reaches nuclear saturation density (applying 
CB parameters as given in Table 1). For simplicity, we assume 
that the CB is expanding with a constant velocity and is moving 
with a constant Lorentz factor. The natural assumption for the 
expansion velocity is the sound speed of the hot blob of matter, 
?;exp ~ Cs ~ Six different cases of CB Lorentz factor 

and mass are investigated (see Table|2j- 

We will first apply a simple model of the CB’s internal en¬ 
ergy obeying a simple equation of state. Using this we calculate 
the evolution of the density, magnetic field and temperature as 
the CB moves away from the origin. In a second step we extend 
our model approach applying an energy equation where pres¬ 
sure degeneracy and neutrino effects are included. As we will 
see, spatial back integration from the transparency radius will 
give strong indication that the CB may be launched close to a 
black hole. The spatial integration back to the source is carried 
out until the CB temperature reaches extreme values, T\ ~ 100 
MeV, unless the CB density reaches nuclear saturation density 
before Tp 

With the CB expanding at a constant expansion velocity 
equal to the sound speed of the matter Uexp = c/ v^, we can use 
certain estimates about the CB at the distance of transparency 
to derive an interrelation between radius and distance from ori¬ 
gin the CB has traveled. As radius of the CB at the distance of 


Table 2. The different cases of CBs explored in this work. 
Note that cases 6, 8 and 9 (Fcb = 1000 and Ncb = 10®^, 
Fcb = 2000 and Nqb = 10®°, Fcb = 2000 and Nqb = 10®^ 
respectively) are not consistent with our assumptions as they 
reaches nuclear saturation density at a distance beyond the light 
cylinder, and have therefore been left out in this paper. 



Fob 

iVcB 

Case 1 

1.0 X 10" 

10« 

Case 2 

1.0 X 10" 

10®° 

Case 3 

1.0 X 10" 

10®® 

Case 4 

1.0 X 10® 

10« 

Case 5 

1.0 X 10® 

10®° 

Case 7 

2.0 X 10® 

10« 


Table 3. The CB radius at the point where it approaches density 
reaches nuclear saturation density {left), at the light cylinder 
Rbc = 1-5 X 10^ cm {middle) and at the distance where the 
CB becomes transparent to radiation {right). 



^CB,nuc [c?7l] 

7?CB,Lc[cm] 

-RcB, trans [c?77.] 

Case 1 

2.5 X 10® 

1.2 X 10® 

1.0 X 10®^ 

Case 2 

5.6 X 10® 

1.2 X 10® 

1.0 X 10®® 

Case 3 

12.0 X 10® 

1.2 X 10® 

1.0 X 10®"* 

Case 4 

2.5 X 10® 

1.2 X lO'* 

1.0 X 10®" 

Case 5 

5.6 X 10® 

1.2 X lO'* 

1.0 X 10®® 

Case 7 

2.5 X 10® 

5.8 X 10® 

1.0 X 10®" 


transparency we apply the estimate by Dado et al. (2002), 

-Rtrans 10^° (g X 10®°) 

To reach this radius, the CB has traveled a period of time 

^ _ -Rtrans f. ,7,7 NcB 

itran^.CB- ~ 577 j sec. (2) 

in the CB rest frame. As it travels essentially with the speed 
of light, at the time when it becomes transparent, the CB has 
traveled a distance 

Tltrans — FcB C ftrans ( 3 ) 

from its origin where it was ejected. This gives a linear scaling 
factor 

^ — T^trans/ Tltrans — 

The radius of a CB is then simply expressed as Rqb = 
D/TcB^/i- The radius at nuclear saturation density and at the 
light cylinder {D\c ~ 1.5 x 10^ cm) are listed in Tabled as 
well as the radius at which the CB becomes transparent to the 
enclosed radiation for the different cases of CB masses and 
Lorentz factors. 
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Fig. 1. Density vs distance from the origin for the CB. The 
backward integration is stopped when the CB density reaches 
nuclear saturation density. 



Fig. 2. Magnetic held vs distance from the origin for the CB. 
The backward integration is stopped when the CB density 
reaches nuclear saturation density. 


Figure ^ shows density as a function of distance travelled 
by the CB. Depending on the CB parameters, the densities at 
the light cylinder range from 10^° — lO^^g/cm^. 

If we estimate the CB magnetic held strength from the 
equipartition assumption, va — Vs = cj with the Alfven 
speed VA and 


B = VA\f^P 



(5) 


3.1. Simple energy equation 

We continue our simple estimates by assuming energy conser¬ 
vation in the CB 


Brad “F F/th T Bmag — Bi^t- (6) 

The radiation energy is written as, 

F'rad = —7ri?Q3 , (7) 

the magnetic energy is, 

Eraag = ^mMcBC^ ( 8 ) 

and the gas thermal energy is, 

Eth = SNcBkT, (9) 

where a = 7.5657 x 10“^^ erg cm“^ K~^ is the radiation con¬ 
stant, k = 1.3807 X 10“^®ergK“^ is Boltzmann’s constant, 
T is the CB temperature, Cm is a parameter that allow us to 
write the magnetic energy in terms of the CB rest mass energy. 
This parameter is hxed by imposing energy equipartition at the 
specihed CB origin. Note that the magnetic energy is constant, 
as we assume that the magnetic held is not dissipated in recon¬ 
nection events and it is not expelled from the CB. The gravita¬ 
tional energy is always negligible compared to the other energy 
channels. The total energy is then given as; 

Etot = 3 X Emag = SCmMcBC^- (10) 

We can now write the energy equation as: 

4/3TrR^aT^ + 3NkT = 2enrMc^ (11) 

This equation is solved to hnd T as a function of D. 

In what follows, we explore two scenarios; (i) the hrst one 
is indicative of the close proximity of a compact source and as 
such it corresponds to the case where the energy equation is 
integrated assuming equipartition at nuclear saturation densi¬ 
ties; (ii) the second rehect scenarios where the CB originates 
from the coronal region of compact stars of of their associated 
accretion disks. Specihcs below. 


3.1.1. Equipartition at nuclear densities: source 
origin 

To set equipartition at nuclear saturation density, Cm has to be 
0.5 for all cases. By rearranging Eq. (O it can be seen that it 
becomes a function of T and p: 

aT'^ + 3pkT/mH = 2pc^em ( 12 ) 


we may compute the CB magnetic held versus the distance 
from origin (Fig.|2j- Close to the hypothetical CB origin where 
we reach nuclear saturation density, the magnetic held strength 
approaches values up to B ~ 10^®G (see also SD- At a dis¬ 
tance from the origin of the order of the light cylinder radius, 
the held strength is between 4.8 x 10^® and 1.4 x 10^® Gauss 
depending on the choice of Lorentz factor and mass of the CB. 
The densities ranges from 10® to lO^^g/cm®. 


The temperature therefore depends on the density only and im¬ 
plies a temperature of about 10^® K at nuclear saturation den¬ 
sity for all cases. The temperature as a function of distance trav¬ 
elled is shown in Fig. Awhile Fig .|4] shows the energy compo¬ 
nents for case 4. The radiation energy is dominant everywhere 
except at nuclear saturation density where there is equiparti¬ 
tion. We note that nuclear saturation density would be reach at 
distances larger than 10® cm for most cases. It is unrealistic to 
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find object with such high densities much larger than 10® cm. 
Also, CBs with such densities need a magnetic field B = 10^® 
G, which is unrealisticly high. We can therefore rule out CBs 
formed with nuclear saturation density. 

It should be noted that the solutions are not very sensitive to 
the choice of total internal energy^. As an example, increasing 
the total internal energy by an order of magnitude we find a 
temperature at the light cylinder, D\c, for case 3 to be 9 x 10^^ 
K, compared to 5 x 10^^ K in our initial calculation. 



Fig. 3. Temperature vs distance traveled by the CB assuming 
equipartition at nuclear saturation density. 



Fig. 4. The different energy components for case 4 using the 
simple energy equation and assuming equipartition at nuclear 
densities. 


^ The expansion energy of the CB is of the same order as the mag¬ 
netic energy. For simplicity we have included it in the expression for 
the total energy. 


3.1.2. Equipartition at the light cylinder: coronal 
origin 

Figure |3 show the evolution of the CB conditions when 
equipartition at a distance of about a light cylinder radius is as¬ 
sumed, and in Fig.|^the energy components for case 4. Table|4] 
shows the temperature found at the light cylinder and the value 
for the energy equipartition parameter for different kine¬ 
matic parameters of the CB (see Tabl^. Note that Cm is now 
determined by the condition that we have equipartition at the 
CB origin (i.e. at a light cylinder distance). 

Table 4. Coronal CB origin. The values for the temperature and 
Cm found for the simple energy equation assuming equiparti¬ 
tion between the magnetic, gas thermal and radiation energy at 
the light cylinder. 


Case 

r[K] 

Em 

1 

5.9 X 10^° 

0.0165 

2 

1.3 X 10“ 

0.035 

3 

2.7 X 10“ 

0.075 

4 

6.1 X 10“ 

0.175 

5 

1.2 X 10“ 

0.325 

7 

1.2 X 10“ 

0.325 



Fig. 5. Temperature vs distance traveled by the CB assuming 
equipartition at the light cylinder. 


By assuming an equal number of electrons rig and baryons 
in the CB the Fermi temperature can be computed. If we take 
the light cylinder distance as a typical length unit at the CB 
origin and compute the Fermi temperature at this location, we 
see that the CB temperature is much smaller than the electron 
Fermi temperature (Tablel^l. We therefore have to improve our 
approach considering also electron degeneracy and neutrino ef¬ 
fects. In the next section we will explore a more appropriate en¬ 
ergy equation where neutrino effects are added. However, we 
will still apply a kinematic approach assuming a linear expan¬ 
sion of the CB. 
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Fig. 6. The different energy components for case 4 using the 
simple energy equation and assuming equipartition at the light 
cylinder. 


Table 5. The values for the magnetic field, temperature and 
density found at the light cylinder for the simple energy equa¬ 
tion. The temperature and magnetic field at nuclear saturation 
density is 1.2 x 10^^ K and 9.3 x 10^^ G respectively for all 
cases, as both the magnetic field strength and temperature is 
dependent on the density only. 


Case B [G] T [K] Tf [K] p [gem"®] 


1 

00 

X 

10 ^® 

5.9 

X 

IQi® 

5.6 

X 

10 “ 

6.1 

X 

10 ® 

2 

1.5 

X 

10 ^® 

1.3 

X 

10 “ 

2.6 

X 

10 “ 

6.1 

X 

10 “ 

3 

4.8 

X 

IQi® 

2.7 

X 

10 “ 

1.2 

X 

10“ 

6.1 

X 

10“ 

4 

1.5 

X 

IQi" 

6.0 

X 

10“ 

5.6 

X 

10“ 

6.1 

X 

10“ 

5 

4.8 

X 

10^'^ 

1.2 

X 

10“ 

2.6 

X 

10“ 

6.1 

X 

10“ 

7 

4.3 

X 

IQi" 

1.2 

X 

10“ 

2.3 

X 

10“ 

4.9 

X 

10“ 


3.2. Coronal CB origin: neutrino effects 


In this section we improve our approach by taking into account 
degeneracy pressure and neutrino cooling (e.g. Popham et al. 
1999). The energy equation becomes. 


Etot + Eiy{t) — Etii + iilrad + T'deg + E^ 

where 

/ \ 1/3 

Edes = SKMcb , 

is the degeneracy energy, 

^ 11 aT^ , ^ 

-^rad = - Mcb, 

4 p 

is the radiation energy and 


the gas thermal energy where 


Xnuc = 30.97 


lOiOg/cm^y ylO^^Ky 


X exp I — 6.096 x 


10i°K 


T 


(13) 

(14) 

(15) 

(16) 


(17) 


gives Xnuc < 1, and Xnuc = 1 elsewhere. In the equation 
above, K = (27r/ic/3)(3/87rm„)^/^ = 1.24 x 10^®, m„ is the 
nucleon mass, R is the gas constant, a is the radiation constant 
and = 2 is the mass per electron. Inserting Eqs. (ID, (HD, 
( I15> and ( I16> into ( I13t gives in terms of D 


Etot + E^{t) = 3 KMcb 


McB3®/2\^/®rcB 


-F -RTMcb 


1 + 3X„ 


47r/ie 

11 




-aT „ 
4 r3 


D 

AttD^ 


(18) 


CB 


35/2 


Note that we apply the same total internal energy (Etot) 
of the CB as in Sec. 3.3, thus the same energy parameter e^a- 
Close to the origin, however, we add an energy component due 
to neutrino effects (emissivity and cooling) and is denoted by 
Efft) inEq.ini 

Two types of neutrino losses may occur, i.e. neutrino emis¬ 
sion due to pair annihilation and neutrino losses due to the cap¬ 
ture of pairs on nuclei. Their contribution to the energy budget 
is computed from Eq. 3.8 and Eq. 3.9 in Popham et al. (1999): 

9 


= 5.0 X 10 


33 


T 


lOii/C 


ergs cm ^ 


(19) 


( 20 ) 

These expressions are integrated over the time it takes the CB to 
reach conditions for which neutrino cooling is not significant. 
We find the latter cooling method (Eq. I20t to be dominant, so 
we limit ourself to using that. 

To compute the effects due to neutrinos, we must know the 
temperature. However, in turn, we want to use the neutrino ef¬ 
fects to find the temperature. We therefore first solve the energy 
equation [O] without adding neutrino effects. Then we use the 
temperature found to calculate the neutrino emissivity which is 
then added to the total energy in equation HD and this equa¬ 
tion is solved to find the temperature as a function of distance 
travelled. The neutrinos are released in small successive bursts, 
mimicking a continuous emission. 

We emphasize again that because of the neutrinos effects 
and the different energy equation used in sec 3.3 there is no 
assumption on energy equipartition applied in this section. 

The temperature is shown in Fig.0as a function of distance 
travelled. In general, the release of neutrinos are seen as a small 
jump in the temperature curve. We note that cases 1 to 4 reach 
the light cylinder with reasonable temperatures (T < 10^^ K) 
and densities (p < lO^^g/cm^). Cases 5 and 7 are ruled out as 
they reach even more extreme conditions before reaching the 
light cylinder. Eor illustrative purposes in Eig.|Dwe show the 
energy components for case 3. Because of the neutrino effects, 
the radiation energy is now the dominant energy, even close to 
the light cylinder. 


4. Sites and formation mechanisms 

In this section we discuss sites that are best suited to account 
for the CBs conditions at the source derived in our previous 
section. We also explore possible formation scenarios. 
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Fig. 7. The temperature assuming the same total energy as in 
sec. 3.3 with neutrino and degeneracy effects. Case 7 goes back 
to 10^^ K (at D = 2.5 x 10^^ cm), whereas case 1, 3 and 4 goes 
to the light cylinder. 



Fig. 8. The energy components for case 3, including degener¬ 
acy and neutrino effects and assuming that case 3 originate at 
the light cylinder. The neutrino release is visible in the radiation 
energy and gas thermal energy at Z? = 3 x 10^ cm. Another 
not visible release occurs at D = 3 x 10® cm. 

4 . 1. “Standard” and hyperaccreting disks 

In annendixlAlwe summarize properties of ’’standard” (cr-disks 
and advection dominated accretion flow disks) and hyperaccre¬ 
tion disks. It is clear that “standard” accretion disks are ruled 
out. Hyperaccreting disks on the other hand are candidates as 
CB sources. They have densities, temperatures and magnetic 
field that are comparable to what we found in the previous sec¬ 
tion. In fact when we consider CBs with hyperaccretion disk 
conditions at the source and perform a forward integration, the 
conditions at Utrans (see Figures lB.llB.3l in appendix|Bji turn 
out to be interestingly similar to those given by Dado et al. 
(2002). The next step then is to look for formation mechanisms 
within the hyperaccreting context. 

The stability of hyperaccreting disks around black holes 
have been recently investigated by Di Matteo et al. (2002) who 


And that the “flows are gravitationally stable under almost all 
conditions of interest”. Exceptions exist for strong accretion 
rates and in the outer part of the disk (see also Narayan et al. 
2001). However as can be seen for Eq. JA.8> - JA.11> in the ap¬ 
pendix these extreme cases favor lower densities and tempera¬ 
tures than those expected for CBs. 

The magneto-rotational instability (MRI, Balbus & Hawley 
1991 and Hawley & Balbus 1991) works only for low mag¬ 
netic held strengths and cannot account for the strong magnetic 
fields required at the origin for CBs. Let us also mention the 
accretion-ejection instability (Tagger et al. 1992) as a possible 
formation mechanism. This instability works for intermediate 
magnetic held strengths and will transfer angular momentum to 
Alfven waves toward the corona of the disk. At extreme magne¬ 
tization the accretion-ejection instability is reminiscent of the 
interchange instability (Spruit et al. 1995) but it seems unlikely 
that these can lead to CBs formation since most of the pertur¬ 
bations are carried by Alfven waves. 

It is thus not clear how a CB can form within a hyperac¬ 
cretion disk. There is also the issue of accelerating the CB to 
r > 100 which is also a major challenge. We will return to this 
in Rafter we discuss other possible formation sites. 

4.2. Neutron tori 

The thick, self-gravitating, neutron tori around 2-3 Mq black 
holes are known to be affected by a runaway instability on time 
scales below the evolutionary time scale of GRBs (Nishida & 
Eriguchi 1996) and we therefore exclude them as source for 
CBs. Simulations of neutron star mergers have also shown that 
about 0.01 Mq of the thick disk of 0.2 Mq around a 1.5 to 3.1 
Mq final central mass distribution becomes gravitationally un¬ 
bound (Rosswog et al. 1999). However, in difference to the hy¬ 
peraccreting disk model, this unbound mass stays rather cold 
(10®iC) and do not constitute a formation site for CBs. 

4.3. Accretion disk corona 

Another possibility is CB formation in the disk corona, for ex¬ 
ample as a huge magnetic flare which ejects a large part of the 
accretion disk corona into a bullet of high velocity. A CB of 
such a size would have a density of about 2.4 x 10® gcm“® 
which is, for comparison, in the range of white dwarf den¬ 
sities. The maximum initial size of the CB we expect not to 
exceed Rlc — 10^ cm. Comparing the CB asymptotic ki¬ 
netic energy to the magnetic energy contained in a volume 
of that size, provides an estimate for the mean magnetic held 
strength of about 10®®G. This corresponds to a magnetic flux 
of ~ 10^® — 10®®Gcm® and is unrealistically high for such 
coronae. 

4.4. Disk-jets and funnei-jets 

Eigure is an illustration of the type of jets that could em¬ 
anate from the vicinity of a compact star. The disk-jet material 
is ejected from the accretion disk while the funnel-jet is ejected 
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(Michel 1969) 


$ 2 ( 72 , 
4Mjetc3 ’ 


( 21 ) 


where <i> = BpV^ is the magnetic flux, ^Ip = c/Dic is the 
angular frequency of the magnetic field and Mjet = TrpVpr^ is 
the mass flow rate within the flux surface. For spherical outflow 
Michel (1969) found that the Lorentz factor at infinity scales as 


roo = CTi/3. (22) 

Fendt & Ouyed (2004) finds a modified Michel scaling in 
the case of a non-spherical magnetic field distribution. In this 
case they find a linear relation between cr and Foo. If the field 
distribution is <I>(r; Af) ~ they find that Foo = 10“^/^cr, 

and if fl)(r; 4') ~ they find 


Fig. 9. Illustration of Funnel-jet and Disk-jet. The funnel-jet is 
launched from a region close to the compact star. The disk-jet 
is launched from the accretion disk. 

from the innermost parts of the disk at the interface with the 
compact star. 

Recent general relativistic magneto-hydrodynamic simula¬ 
tions by De Villiers et al. (2005) of a black hole and an ini¬ 
tial torus seeded with a weak poloidal magnetic field show that 
a funnel jet with Fee ^ 50 is formed. Instabilities do occur 
in funnel-jets, however, the induced instabilities have densities 
much lower than the CB values found in the previous section. 
Funnel-jets can therefore be ruled out as a possible formation 
site for CBs. 

A disk-jet becomes cylindrically collimated on a length 
scale of the order of 1-2 light cylinder distances (Fendt & 
Memola 2001). Knot generating instabilities reminiscent of 
CBs are known to occur as jets collimate (Ouyed et al. 1997). 
This is a possible formation mechanism for CBs. What remains 
is to show how they can be accelerated to high F. Specifics be¬ 
low. 

5. CB acceleration to ultra-relativistic velocity 

Having isolated, or more precisely favored, jets from hyper¬ 
accretion disks as plausible formation sites for CBs we now 
discuss acceleration mechanisms with which CBs can reach 
Lorentz factors in the thousands. 

Assuming that the CB is accelerated by converting the in¬ 
ternal magnetic energy to kinetic energy, we can find an es¬ 
timate for the magnetic field needed to explain such Lorentz 
factors. Using typical CB radii close to the source and parame¬ 
ters from Table [Owe find that the magnetic field must be of the 
order of 10^® — 10^® G. This is unrealistically high. Either our 
approach is too simple, or a different acceleration mechanism 
must be at work. 

5.1. MHD Acceleration: CB speed 

The ability of the magnetic field to accelerate particles to 
high Lorentz factors is given by the magnetization parameter 


Foo = 10-i/V , (23) 

in which case hyperaccreting disks with ejection rates of the 
order 1Q~^Mq/s and magnetic field of the order 10^'^G can 
lead to jets with a Lorentz factor Foo ~ 1875. 

5.2. MHD instability: CB mass 

To a first order, instabilities related to Alfven crossing time can 
develop on timescales 

fi„s = tA = ^, (24) 

VA 

where is the radius of the disk-jet. For li?ic < 77jet ^ 
10i?ic, we arrive at tins ~ 1 — 10 ms which would imply the 
plausible formation of blob of matter as massive as Mjns = 
tinsMjet ~ 10“® — 1O“^M0. This can be compared to the 
typical CB mass of the order Mcb = 10“^Mq. 

As we have shown above, first forming the CB in the disk 
and then accelerating it will require unrealistic magnetic fields 
of the order 10^® G. However, first accelerating the wind to the 
light cylinder and then forming the CB through an instability 
beyond the light cylinder requires much smaller magnetic field 
strength (< 10^"^ G). This is a possible mechanism for forming 
and accelerating CBs. 

6. Conclusion 

Assuming that CBs move and expand with a constant veloc¬ 
ity we have estimated the CB conditions as close as possi¬ 
ble to their origin. CBs require extremely high internal mag¬ 
netic fields when they are formed with field strength exceed¬ 
ing ~ 10^® G. The temperature was found to be of the order 
of 10^^ — 10^® K. The physical parameters of the CBs at the 
origin are, within an order of magnitude estimates, indicative 
of hyperaccreting disks. However, if formed in the accretion 
disk we find it challenging to accelerate the CBs to the high 
Lorentz factors. The coronal origin is ruled out because of the 
unrealistically high coronal magnetic flux necessary to form the 
CBs. Our results instead hint to a jet origin of CBs. The radius 
(< Die) and mass flow ( s) in a jet from a hyperac¬ 

creting can account for the CB mass and density. Furthermore, 
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this outflow can be accelerated to F ~ 2000 by MHD pro¬ 
cesses (Fendt & Ouyed 2004). Any instability in this outflow 
beyond the light cylinder could lead to CB formation. We thus 
suggest that CBs form as instabilities in ultra-relativistic jets 
emanating from the surface of hyperaccretion disks. The tight 
link between SNe and the CB model for GRB requires that all 
(or almost all) core collapse SNe will produce CBs. Our work, 
within its limitations, implies that hyperaccretion disks must be 
a common occurrence in core collapse SNe to accommodate 
the CB model - a notion which remains to be confirmed. 
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Appendix A: Accretion disks 

A.1. “Standard” accretion disks 


assuming the default parameters. For the equipartition mag¬ 
netic held one gets 


B[G] = lO*"! 
The temperature is 

T[K] = 1.3 X 10 


/ M 


\M(, 


- 1/2 


fe) 


-3/4 


(A.2) 


0.001 


-1/4 




-3/4 


(A.3) 

Advection Dominated Accretion Flow (ADAF) disks have 
density (Narayan & Yi 1995): 


p[gcm = 6.5 X 10’ 

- 1-1 


-3 I 


f—y 

V AfEdd / 


0.001 

-1 




Cl Cg 
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where ci and C 3 are defined in Eq. (2.1) in Narayan & Yi 
(1995). 

The corresponding magnetic held is 


B[G] = 5.5 X 10' 

J- 


- 1/1 


0.001 
M 




) VMeJ 


V A^Edd 


-5/4 


(A.5) 


while the ion temperature of such disks are (Narayan et al. 
1998) 

r 

, 2 rs 


Ti[K] = 2 X 10^2/3 
where (3 is given by 


(A. 6 ) 


= ^ = (1 - 


(A.7) 


A.2. Hyperaccreting disks 

Hyperaccreting disk (Popham et al. 1999) density is: 

-1-3 / \ 4-1 


p[gcm = 1.3 X 10 


12 I 
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M 


M 


MqS- 
-1.7 / \ -2.55 
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(A. 8 ) 


their disk scale height is 


A standard Shakura-Sunyaev disk (Shakura & Sunyaev 1973) 
will have a density 


p[gcm-3] = 7.2 X lO-'^ 


Cky 

0.001 


( 


3rs 


3/2 


X 




(A.l) 


where is a viscosity parameter, rs is Schwarzschild radius 
and MEdd is the Eddington mass. With this density, the radius 
of a CB with mass M = 10^° baryons becomes 3.8 x 10® cm. 


iT[cm] = 1.9 X 10*^ 
while the temperature is 
T[K] = 7.6 X 10 
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The corresponding equipartition magnetic held is of the order 
of 

B[G] ~ 10^^ - 10^®. 


(A.ll) 
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Appendix B: Forward integration 

For completeness and for self-consistency check, here we con¬ 
sider CBs with hyperaccreting disk conditions at the origin and 
perform a forward integration until the CBs reach the distances 
where they become transparent. 

Assuming that the CB radius evolves as before, R = 
c/ (rcB>/3), then the density at a distance corresponding to the 
surface of the hyperaccreting disk {D = 10^ cm) will be too 
high. We will thus make a slight adjustment by rewriting the ra¬ 
dius as i? = D/{Tcb'</^)+x, where a; is a number that ensures 
that the density at the origin does not exceed p = lO^^g/cm^. 
Therefore, x is found by solving the following equation; 


Mqb 


ctxR = --- 


which implies: 


a: [cm] = 


(6/7r)i/3rcB - 712 X 10^ 


GOOOOFcb 


(B.l) 


(B.2) 



Fig. B.l. Magnetic field vs distance from the origin for the CB, 
starting with hyperaccreting disk conditions. The forward in¬ 
tegration is stopped when the CB becomes transparent to its 
enclosed radiation. The initial magnetic field for all cases is 
about 5 X IQi® G. 


This also ensures the correct expansion velocity Uexp = c/v/S- 
Table IB. ll shows the corresponding parameter values for x and 
e„i- We should also note that in this case e„i will be chosen as to 
insure energy equipartition at the disk surface. The temperature 
thus found is used to calculate the neutrino emissivity, which is 
then added to the total energy in Eq.[21to find the new temper¬ 
ature. As before, the neutrinos are released in small successive 
bursts mimicking a continuous release of neutrinos. 

Figs. IB.1IB.3I shows the temperature, magnetic field and 
density as a function of distance. All cases starts with T ~ 
2 X IQii K, p = 10i2g/cm3 and B = 6 x 10^® G at D = 
10® cm. The neutrino effects can be seen as small jumps in the 
temperature curves, but in general the neutrinos do not change 
the overall picture a lot. The neutrino contribution were of the 
same order or smaller than the total energy, and as discussed 




Fig. B.l. Temperature vs distance from the origin for the CB, 
starting with hyperaccreting disk conditions. The forward inte¬ 
gration is stopped when the CB becomes transparent to its en¬ 
closed radiation. The initial temperature for all cases is about 
2.7 X 10“ K. 


Fig. B.3. Density vs distance from the origin for the CB, start¬ 
ing with hyperaccreting disk conditions. The forward inte¬ 
gration is stopped when the CB becomes transparent to its 
enclosed radiation. The initial density for all cases is about 
1 X 10®^g/cm®. 

Table B.l. The parameter x used in the relation between Rcb 
and D, and for the different CB cases when starting with 
disk conditions and integrating forward. 


Case 

X [cm] 

Cm 

1 

15279.2 

0.045 

2 

33584.5 

0.045 

3 

73022.3 

0.045 

4 

15798.9 

0.045 

5 

34104.2 

0.045 

7 

15827.7 

0.045 


before the temperature is not very sensitive to changes in the 
total energy. 
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The temperature at Dtrans is of the same order as for the 
backward integration (T = 10"^ K to T = 10® K), and also 
close to the value given by Dado et al. (2002) of Ttrans — 4 eV. 
The difference between the backward and forward integration 
at large distances is due to the different Cm parameter. For large 
distances, the a:-parameter does not play any role. 

To summarize, the results of the forward integration indi¬ 
cate that CBs formed within hyperaccretion disks could in prin¬ 
ciple provide the necessary conditions at Utrans to account for 
GRB features as claimed in Dado et al. (2002). 


